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Abstract 

In the standard model there are charges with abelian anomaly only (e.g. 
right-handed electron number) which are effectively conserved in the early 
universe until some time shortly before the electroweak scale. A state at 
finite chemical potential of such a charge, possibly arising due to asymme- 
tries produced at the GUT scale, is unstable to the generation of hypercharge 
magnetic field. Quite large magnetic fields (~ 10 22 gauss at T ~ 100 GeV 
with typical inhomogeneity scale ~ ^§r) can be generated. These fields may 
be of cosmological interest, potentially acting as seeds for amplification to 
larger scale magnetic fields through non-linear mechanisms. Previously de- 
rived bounds on exotic B — L violating operators may also be evaded. 
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It is usually assumed that the early Universe at temperatures above the electroweak scale 
and below, say, 10 12 — 10 16 GeV (depending on the model of inflation) consists of an (al- 
most) equilibrium primordial plasma of elementary particles, in which any long-range fields 
are absent. One exception is in the context of the problem of generating galactic magnetic 
fields, which may require the presence of primordial seed magnetic fields which are subse- 
quently amplified by a galactic dynamo mechanism (see, e.g., 0). The creation of long 
range magnetic fields requires that conformal invariance be broken in the coupling of the 
electromagnetic field to gravity 0, and a number of mechanisms based on different ideas 
about this breaking have been proposed to date ||[| . In this letter we argue that there may 
be a relation between the appearance of magnetic fields in the early Universe and two other, 
apparently completely unrelated, phenomena : (i) The smallness of the electron Yukawa 
coupling constant, and (ii) possible lepton asymmetry of the early Universe. 

In short, the logic goes as follows. There are three exact conservation laws in the standard 
electroweak theory. The associated conserved charges can be written as N$ = Li — where 
Li is the lepton number of iih generation and B is the baryon number. The fourth possible 
combination, B + J2iLi is not conserved because of electroweak anomalous processes, which 
are in thermal equilibrium in the range 100 GeV < T < 10 12 GeV ||. Now, if h e = 0, 
where h e is the right electron Yukawa coupling constant, then the electroweak theory on 
the classical level shows up a higher symmetry, associated with the chiral rotation of the 
right electron field. For the small actual value of the Yukawa coupling (h e = 2.94 x 10 -6 
in the MSM) this symmetry has an approximate character. At temperatures higher than 
Tr ~ 80 TeV perturbative processes with right electron chirality flip are slower than the 
expansion rate of the Universe ||, and therefore this symmetry may be considered as an 
exact one on the classical level at T > Tr [[J. (The importance of this symmetry for the 
consideration of the wash-out of the GUT baryon asymmetry by anomalous electroweak B 
and L non-conserving reactions was realized in refs. 0,^|,§|.) Suppose now that an excess 
of right electrons over positrons was created by some means at T > Tr ( e.g. by a GUT 
mechanism for baryogenesis) . Now the right electron number current is violated in the 
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minimal standard model (MSM) as described by the anomaly equation 

where /(/) are the ?7y(l) hypercharge field strengths (and their duals) respectively, g' is 
the associated gauge coupling and Ur = —1 is the hypercharge of the right electron. The 
number of the right electrons Nr therefore changes with the Chern-Simons (CS) number of 
the hypercharge field configuration as ANr = ^y R AN cs with 

J2 



32tt : 



J d 3 xe ijk fijb k , (2) 



where b k is the hypercharge field potential. 

One can now see qualitatively that there is an instability in hot matter with an excess 
of right electrons towards formation of hypercharge fields with CS number as follows. (The 
line of reasoning presented here is similar to the consideration of cold fermionic matter with 
anomalous charges in |J). The energy density "sitting" in right electrons with a chemical 
potential hr is of order UrT 2 , and their number density of order [irT 2 . On the other 
hand this fermionic number can be absorbed by a hypercharge field of order g' 2 kb 2 , with 
energy of order k 2 b 2 , where k is the momentum of the classical hypercharge field and b is its 
amplitude. Therefore, at b > T/g' 2 and k ~ \irT 2 / (g' 2 b 2 ) the gauge field configuration has 
the same fermion number as the initial one, but smaller energy. An instability to generation 
of hypercharge magnetic field, which tends to "eat up" real fermions, results. It is important 
here that at temperatures T > Tr the electroweak symmetry is "restored", and that the U(l) 
hypercharge magnetic field is massless at that time. (No term like myb 2 is generated in any 
order of perturbation theory in abelian gauge theory at high temperature [fL0[| ; the lattice 



study in [[Tl|] confirmed this expectation for SU(2)xU(l) EW theory beyond perturbation 
theory). If the hypercharge magnetic fields survive until the time of the EW phase transition 
(T ~ 100 GeV), they will give rise to ordinary magnetic fields because of electroweak mixing. 
In the rest of this paper we present quantitative estimates of the (hypercharge) magnetic 
fields which may be produced by this effect. 



Let us discuss first the possible origin and the magnitude of the the required right 
electron number asymmetry 5r = e^/s, where s = ^Ti 2 T 3 N e ff is the entropy density with 
N e ff = N), + ^Nf = 106.75 the total effective number of degrees of freedom of the MSM. 
In principle 5r produced by out of equilibrium decay at the GUT scale can be as large as 



the final baryon asymmetry 5 b being that observed since there is no simple general relation 
between the two numbers. In theories like those discussed in [|,|J with L violating processes 
at intermediate scales one has 5b ~ 5r, at least in the case that the L violating processes 
go out of equilibrium before the cr violating ones come into equilibrium. In || the case is 
considered where the L violation continues for just long enough to reduce the final 5b to the 
observed one from an initially larger value fixed by 5r. And, in a simple GUT like SU(5) in 
which the charges iV, = Li — |5 are conserved, we can have 5b = at the electroweak scale 
irrespective of the value of 5r during the time it is effectively conserved. In the rest of this 
letter we will simply assume the existence of a primordial density of e#, with its chemical 
potential as a free input parameter, assumed only small enough to be treated perturbatively. 
We also assume that no hypercharge magnetic fields existed before the right electron excess 
is generated. 

The effect of the anomaly on the gauge field dynamics is given through the term in the 
effective Lagrangian 



which is obtained by integrating out the fermions at finite chemical potential |13| . It simply 
describes how winding the gauge fields to give CS number changes the energy of the system 
because it changes the number of fermions as described by the anomaly equation. Adding 
this term to the ordinary Lagrangian for the gauge fields leads to the equations of motion 



10-2 _ iQ-4 



(for a review see, e.g. |I2||). This is quite consistent with the magnitude of 



5C 




pR^ijkfijbk- 



(3) 



where a 1 = g' 2 /4ir. These are simply Maxwell's equations with the additional term due to 
the anomaly (oc HrH) and the assumption that the total (hypercharge) current is given by 
j = crE, where a is the conductivity of the plasma and E is the (hyper-)electric field. We 
have also dropped the term BE / dt since (as we will see below) the fields always evolve on a 
time-scale which is much longer than <r -1 . In the expanding FRW Universe with scale factor 
a the equations have exactly the same form in conformal time coordinates r = / a~ 1 (t)<it, 
but with the replacements /j,r — > fiRa and <r(oc T) — > oa. The fields E and H are those given 
by their standard definitions in the conformal frame which will be related to the physical 
fields at the appropriate point below. We also have the following kinetic equation for hr : 

ld(/i R a) a' 783 1 - - 

- a -¥T = T¥^ £ ' H ~ VrM ' (6) 

in which the first term describes the change in the chemical potential due to the anomaly 
(// oc E-H), and the second the change due to the perturbative processes which flip electron 



chirality with the rate T R = j^T (M = M p i/1.66^N eff ~ 7.1 x 10 17 GeV). The numerical 
coefficient ^ comes from the relationship between right electron chemical potential and 
right electron number asymmetry (in terms of which the anomaly is expressed) 
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which is obtained from a local thermal equilibrium calculation in the EW theory with three 

fermionic generations and one scalar doublet, with the conserved charges assumed to be 

Ni(= 5is) and e R (= S R s). 

With a Fourier mode decomposition H(x) = f d i kH(k)e~ lk ' x with H(k) = h^i where 

% = 1, 2, el 2 = 1, e*j ■ k = 0, e\ • e.2 = 0, the linear equations @ and (|5]) become 

_ , k 2 Aifi\k\ . 

d T h x H hi ! - L - -h 2 = 

aa a 

k 2 , 4:iu\k\ , n . , 

d T h 2 + —h 2 + — = (7) 
era a 

where \i = f^fiR- The mode 

/i 2 (r, k) = -ih^r, k) = -^(/ii(0, fc) + z/i 2 (0, fe))ezp(A+(r)) (8) 



where 

A±(r) = -— (Jfer T 4 T rfr'ua) (9) 
aa jo 

is an unstable mode which is growing at conformal time r if k < 4/i(r)a(r). It has the 
property E{k) = —(—\k\ + 4fia)H(k). (The other orthogonal mode decays at any \k\.) 

Consider now the approximation in which the chemical potential /i is a constant. The 
growing instability starts to develop at T ~ T g where we define T g to be 

8(i)»'^ = l. (10) 
T a/T T g 

(when the maximally growing mode with k = 2fia has begun growing significantly). A nec- 
essary requirement for the instability to develop is that T g > T R , since if this is not satisfied 
the second term in (^) will rapidly reduce towards zero. Translated into a minimum value 
for 5 = £ this requires 6 > 5 cr i t = 1CT 6 (using o « 68T WW). For 5 < 6 cr u no non-trivial 
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dynamics result from the presence of such a chemical potential since the unstable modes are 
frozen on the relevant time-scale. If 5 > 5 cr i t the evolution of the instability for T <T g will 
be given by the simple growth factor above, until the time at which the growth becomes 
significant enough that the first term in (||) is important. To estimate when this is and what 
the amplitude of the field is at that time it is sufficient to calculate the CS number as a 
function of time. It is given (per co-moving volume) by 

g' 2 

n cs(r) = -^-^ < e ijk f jk {r)bi(T) > 
62,71 

neglecting all but the growing mode. We have also taken < bi(k, r)b*(l, r) > | T=0 = 
5 3 (k — l)5ij < b 2 (k) >o, assuming translational and rotational invariance of the initial per- 
turbations, and assumed that the perturbations are thermal in origin, with the appropriate 
normalizations, < b 2 (k) >o= /(A:) where f(k) = (e^o — 1) _1 is the bosonic distribution 

function and Tq the temperature at which we define clq = 1. Defining e by \y 2 R n cs = eAe R a 3 , 
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where Ae# = ^^T 3 i.e. the difference between the right electron density in the initial 



state and the fin = state, the linear approximation breaks down when e ~ 1. Evaluating 
the integral in (|TT| ) we find esi2x l0~ 6 5-^e a where a = T g /T with T g as in (|T0|). Thus 
for a few expansion times after the mode starts growing at temperature T g we have e < 1 
and the linear approximation is valid. The corresponding physical magnetic field H p h y can 
be estimated by putting |n cs | m -^-jkb 2 = |eAe^a 3 and using kb = a 2 H phy} where k ~ 2/xa 
(i.e. assuming the maximal growing mode to dominate). Putting in the numbers this gives 



a physical magnetic field of strength H P h y ~ 2 x 10 2 \J tS^r-T 2 at a physical length scale 
kpL = ~ Wr- For 5 ~ 5 x 10~ 6 we are in the linear regime until Tr ~ 80 TeV and 



at that point therefore have a magnetic field of strength H p f iy ~ 6 x 10 26 gauss (1 GeV 2 



v phy — \ a J 25T- 

• i „ i-:„ c„u „£ „j. ii. l 

l phy 

= 1.95 x 10 20 gauss) at a length scale of ~ 10 5 /T (compared to a horizon scale of ~ 10 13 /T). 

How do the fields evolve for T < Tr? In the case that the linear ( i.e. constant chemical 
potential) approximation is good, one expects that the growth will rapidly turn into decay 
as fi is damped. Within a few expansion times the growth will be undone as the maximally 
growing mode k ~ 2fia now decays with exponent —(k 2 /a)ar. What about the case when 
this linear approximation breaks down? To treat this case we must analyse the full non-linear 
set of equations ([§)-(§[)■ We have done this numerically with the simplifying assumption 
that the distribution of right electron number is homogeneous in space. Then the two linear 
equations (§]) and (H) can be solved exactly for any time dependent /z, the solution inserted 
in (H), and the averaging over thermal initial conditions performed. The resulting equation 
is 

Bun <V 783 1 r°° 

-p = (_)*_-__ / dkk 2 [(k - ^a)exp(2\ + (r)) - (k + 4/xa)exp(2A_(r))], (12) 
or ix 88 6Z7i A aa z l Jo 

a length scale of 10 5 /Ta length scale of 10 5 /Twith A±(r) given by @. Our results show 
that the chemical potential, typical physical momentum of the magnetic field configuration 
and the magnetic field energy scale as 

L>>r k p h y H h /T\ 2 

T T T 4 \T g ) v ; 

in the range T g > T > Tr. This behaviour can be easily understood qualitatively as follows. 
As the instability develops, the linear approximation breaks down and [i starts significantly 



decreasing. This shifts the growth of modes to longer wavelengths. This procedure continues, 
growth of any mode eventually turning itself off and increasing the growth coefficient of 
modes at larger scales. The minimum value of \x which can be reached at any given time 
r (and, correspondingly, the maximum physical scale for the sourced fields) is simply that 
given by ([IT]), solved for \i with T g replaced by the temperature T(r) i.e it is just the minimal 
chemical potential required to drive a growing mode at that time in the linear approximation. 
The parametric dependence on the temperature observed follows from the fact that the 
chemical potential (and maximally growing mode) trace these values. The dependence of 
the magnetic field energy follows from the expression we derived in the previous paragraph 
by setting the CS number of the configuration to cancel the total fermion number, but now 
taking the appropriate scaling for k p h y itself. 

Evolving the system forward from Tr to the electroweak scale T ew we see the damping 
of the fields for 5 in the linear regime anticipated above, as the perturbative processes erase 
the chemical potential driving the growth. As 5 increases, however, this damping becomes 
less efficient, and for 5 > 2 x 1CT 4 we find that the damping has not set in at all by the 
electroweak scale. The reason for this behaviour is also simple. For a mode which evolves 
in the linear regime, the growth and decay exponents are effectively the same for the modes 
which grow significantly. Once we enter the non-linear regime this is no longer true, since 
the maximally growing mode carries the integrated effect of growth until any given time i.e. 
it has grown with exponent k J ^-dr which is much greater than — r. Put another way, 
the mode has been able to grow on a scale significantly larger than the diffusion length for 
magnetic field at the relevant time, and it takes some time after the end of the growth for 
the latter scale to catch up and undo the effect of the instability. For 5 in this region we also 
see that the typical scale of magnetic fields k and value of the chemical potential at T ew do 
not depend on initial asymmetry, —^m. ~ ^ ~ The amplitude of magnetic field scales 
as H oc y/6 and e.g. for 5 = 1CT 2 we find H ~ 4 x 10 22 gauss. This is as we would expect 
from the discussion of the scaling above. 

What is the ultimate fate of these magnetic fields? Unless some other effect comes into 
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play in the dynamics, the fields will decay. One such effect is turbulence. With the full 
set of MHD equations (which include the velocity of the fluid which we have neglected) 
there is a transition to a turbulent regime when the magnetic Reynold's number R = aLv 



is large [[TIJ. The reason we have evolved the equations to the electroweak scale is that, if 
the electroweak phase transition is of first order, it serves as a source of turbulence [ITS . 
Since we have here a ~ 10 2 T and magnetic fields which begin to grow on length scales L 
up to ~ 2 x 10 6 /T we expect to enter the turbulent regime if there are bulk velocities of 
greater than ~ 5 x 1CT 9 , which are certainly larger that the expected velocity of the bubble 
walls. A recent study of this phenomenon |l6j suggests that the effect of this turbulence is 



to transfer the magnetic energy to larger length scales, thus evading the Silk argument |]T7 
If true, the fields generated by the mechanism under discussion may play the role of the 
seed galactic magnetic fields. Note that the seed fields we obtain at the electroweak scale 
with the mechanism we have discussed (~ 10 22 gauss) are much larger than those generated 
at bubbles walls (~ 1CT 2 gauss) which were suggested as seeds for amplification through 
turbulence in [JTJJ. It is also worth mentioning the particular structure of the magnetic fields 
appearing because of the abelian anomaly. The CS wave (^) has a non-zero value of H -d x H 
and thus breaks parity. Could it be that the rotation of galaxies are related to this? Study 
of the entire set of MHD equations with the additional anomalous terms discussed in this 
letter will be required to address this question. 

Finally let us mention that the processes we have considered also affect the bounds on 
the strength of exotic interactions with B — L violation derived from the requirement that 
GUT baryon asymmetry is not erased by sphalerons |5]|| (which is important if no baryon 
asymmetry is created at the electroweak scale). If the right electron asymmetry produced 
at the GUT scale is small enough (5 < 1CT 6 ), then the bounds are obviously not affected 
since abelian anomaly does not play any role. If, on the other hand, 5 is large enough that 
significant CS number survives remains in the condensate until the electroweak scale, any 
bound on the strength of exotic interactions can be evaded. Irrespective of the effect of any 
B — L violation until that point the remaining CS number will be converted into quarks and 



leptons carrying net baryon number at the electroweak phase transition. The final baryon 
asymmetry will depend on the initial value of 5 and the exact strength of the phase transition 
(which will determine how the B violating processes turn off). Conversely, given detailed 
knowledge of the phase transition, it will be possible to place an upper bound on the initial 
value of 5 in the very early universe, and on the strength of the magnetic fields resulting at 
the electroweak phase transition. 

We are grateful to K. Enkvist, M. Giovannini, A. Kusenko and L. McLerran for interesting 
discussions. 
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